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摘 要 :针对 无 人 机 集群 截获 的 信号 样本 难以 直接 融合 分 析 , 以 及 训练 样本 较 少 且 工作 模式 样本 不 平衡 条 件 下 多 功能 雷达 
(MFR) 工 作 模 式 识别 精度 低 的 问题 ,提出 了 一 种 基于 自 组 织 映 射 (SOMD) 聚 类 平滑 图 信号 生成 的 MFR 工作 模式 识别 方法 。 
首先 ,利用 分 布 式 SOM 算法 对 截获 的 信号 样本 集 进行 聚 类 ,提取 样本 之 间 的 相似 性 ;然后 ,依据 聚 类 结果 将 信号 样本 集 以 
平滑 图 信号 的 方式 表征 ,建立 同一 工作 模式 下 信号 样本 的 关联 ;最 后 ,采用 图 注意 力 网 络 对 上 述 图 信号 进行 图 节点 数据 融 
合 与 分 类 ,完成 MFR 工作 模式 识别 。 实 验 结果 表明 ,在 工作 模式 样本 不 平衡 度 约 为 10:1 ,每 种 类 别 训练 样本 数 为 25 时 ,该 
方法 的 识别 准确 率 和 F1 指数 相对 现 有 方法 分 别提 高 了 22. 896,22. 34%, 且 能 适用 于 存在 一 定 噪声 干扰 的 情况 。 
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Abstract: UAV swarms are widely used in radar signal interception due to their advantages of wide sensing range and rapid 
: information sharing. Aiming at the problem that the signal samples intercepted by UAV cluster are difficult to be fused and 
lyzed directly. and the recognition accuracy of multi-function radar (MFR) working mode is low under the condition of 
few training samples and unbalanced working mode samples. an MFR working mode recognition method based on smooth 
2 mph signal generated by self-organizing map (SOM) clustering is proposed. Firstly. the intercepted signal samples are 


stered by using distributed SOM algorithm to extract the similarity between samples; Then, according to the clustering 


ulti, the signal sample set is characterized by smooth graph signal. and the correlation of signal samples under the same 
working mode is established; Finally. the graph attention network is used to fuse and classify the graph node data of the a- 
bove graph signals to complete the MFR working pattern recognition. The experimental results show that, when the imbal- 
ance of working mode samples is about 10:1 and the number of training samples in each class is 25, the recognition accuracy 
and F1 measure of this method are improved by 22. 8% and 22. 34% respectively compared with the existing methods. and 
can be applied to the case of noise interference. 
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雷达 工作 模式 识别 是 电子 侦察 的 核心 内 容 , 其 通 ”的 工作 模式 ,为 制定 电子 对 抗 策略 提供 情报 基础 。 
过 对 截获 的 信号 样本 进行 分 析 处 理 以 判定 目标 雷达 ”多 功能 雷达 (multi-function radar. 简称 MFR) 因 具 
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有 多 工作 模式 、 智 能 程度 高 可靠 性 强 等 优点 被 广泛 
应 用 于 电子 战 ,给 雷达 信号 侦察 带 来 巨大 挑战 。 
MFR 信号 被 截获 概率 较 低 ,因此 在 短 时 间 内 可 截获 
的 信号 样本 较 少 ; 若 不 对 MFR 施加 一 定 对 抗 压力 ， 
其 某 些 工作 模式 可 能 被 隐藏 沾 。 将 无 人 机 集群 应 用 
于 MFR 侦察 可 有 效应 对 上 述 问 题 。 其 不 仅 能 通过 
对 战场 的 全 覆盖 与 信息 交互 得 到 较为 完备 的 MFR 
信和 号 样本 集中 ,还 可 实施 模拟 攻击 对 MFR 施 压 ,使 
其 暴露 所 有 工作 模式 , 且 应 用 无 人 机 集群 是 电子 对 抗 
REKKAA”, 

在 复杂 电磁 环境 下 ,截获 的 信号 样本 通常 存在 脉 
WEAK ,测量 误差 较 大 等 缺陷 '* ,通过 数据 融合 能 够 
减轻 其 对 识别 结果 的 影响 。 然 而 ,由 于 各 无 人 机 截获 
的 MFR 信号 样本 通常 来 自 多 个 MFR 或 不 同 工 作 模 
式 5 难 以 直接 对 其 进行 融合 处 理 。 因 此 ,如 何 建立 无 
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号 样本 通常 具有 一 定 相似 性 的 特性 ,生成 关于 信号 样 
本 集 的 平滑 图 信号 ,图 信号 中 存在 连接 的 信号 样本 通 
常 具 有 和 较 强 的 相似 性 , 以便 对 信和 号 样本 进行 融合 处 
理 ; 利 用 图 注意 力 网 络 (graph attention networks, 简 
称 GAT) 对 上 述 平滑 图 信号 进行 一 阶 邻 域 数据 融合 
与 节点 分 类 ,实现 工作 模式 样本 不 平衡 且 训 练 样本 较 
少 情 况 下 的 MFR 工作 模式 识别 。 


1 场景 描述 

考虑 的 雷达 对 象 为 “水 星 ”MFR, 它 是 一 部 真实 
装备 的 防空 雷达 ,Visnevski 虽 未 公开 该 雷达 的 信和 号 
参数 .工作 模式 转换 规律 等 信息 ,但 给 出 了 行为 层面 
的 特征 。“ 水 星 ”MFR 采用 符号 语言 建 模 一 种 多 层级 


人 机 集群 截获 的 MER 信号 样本 之 间 的 融合 依据 是 
必须 考虑 的 问题 。 

TAAK MER 工作 模式 识别 方法 主要 包括 基于 
[ze 尔 科 夫 模型 (hidden markov model, 简称 
HMMD 的 方法 、 基 于 句法 模型 的 方法 、 基 于 信号 特征 
提取 的 方法 。 基 于 HMM 的 方法 5 将 雷达 系统 视 
X t oos 自动 机 (finite state automaton ,简称 
FSA) ,利用 HMM 对 FSA 建 模 ,实现 MER 工作 模 
式 讽 别 。 该 方法 在 保留 观测 信和 号 最 大 信息 量 的 同时 
具有 很 好 的 灵活 性 ,但 难以 对 状态 数 繁多 、 转 换 复杂 
的 震 达 系统 进行 准确 表征 。 基 于 句法 模型 的 方 
VO didt MER 工作 模式 转移 的 文法 类 型 ,构造 相 
应 苦 自 动机 , 实现 工作 模式 识别 。 该 方法 理论 性 较 
强 ? 引 相 较 于 基于 HMM 的 方法 具有 更 好 的 雷达 信 
号 恒 模 能 力 ,但 它 十 分 依赖 雷达 的 先 验 知识 ,本 质 上 
是 对 情报 库 中 先 验 知识 的 符号 化 表征 。 基 于 信号 特 
征 提 取 的 方法 nm" 依据 不 同 工 作 模式 下 信号 特征 的 
差异 性 ,构建 分 类 器 以 实现 MER 工作 模式 识别 5 。 
该 方法 对 先 验 知识 的 需求 较 低 ,识别 准确 率 通常 较 
高 , 且 对 含 噪 数据 具有 较 强 的 适应 性 ,但 训练 复杂 分 
类 模型 所 需 的 大 量 训练 样本 不 易 获 得 。 此 外 ,由 于 


信号 结构 , 即 雷达 脉冲 .雷达 字 、 雷 达 短 语 。 雷 达 短 语 
由 一 定数 量 的 雷达 字 串 联 而 成 ,雷达 字 由 一 定数 量 的 
雷达 脉冲 串联 而 成 52 。 此 外 ,雷达 短语 之 上 还 有 雷 
达 句 子 层级 ,但 不 在 本 研究 范围 内 。 

“水 星 ”MFR 采用 雷达 短语 表达 工作 模式 ,有 搜 
索 . 捕 获 、 非 自 适应 跟踪 、 距 离 分 辨 .跟踪 保持 5 种 工 
作 模 式 。 由 于 该 雷达 存在 少量 复 用 的 雷达 短语 ,而 本 
方法 未 考虑 该 问题 ,因此 对 部 分 雷达 短语 作 了 归并 处 
理 , 处 理 后 的 “水 星 ”MFR 工作 模式 及 其 对 应 的 雷达 
短语 如 表 1 所 示 。 

对 抗 场景 如 图 1 所 示 , 其 展示 了 MFR 工作 模式 
识别 的 一 般 流 程 。 首 先 ,利用 无 人 侦察 机 截获 MFR 
的 信号 脉冲 流 ; 然 后 ,经 过 参数 测量 、 信 和 号 分 选 、 雷 达 
字 提 取 等 处 理 提 取出 雷达 字 序 列 ; 最 后 ,从 雷达 字 序 
列 中 切割 出 一 个 个 雷达 短语 ,对 其 进行 分 类 ,以 实现 
MFR 工作 模式 识别 。 本 研究 的 处 理 对 象 是 雷达 短 
语 , 而 参数 测量 ,信号 分 选 .雷达 字 提 取 等 算法 不 在 本 
研究 范畴 内 ,但 考虑 了 因 上 述 算法 存在 误差 引起 的 信 
号 样本 缺陷 ,具体 而 言 ,在 仿真 中 考虑 了 存在 测量 误 
差 . 脉 冲 丢失 和 干扰 脉冲 的 情况 。 

对 无 人 机 集群 截获 的 雷达 短语 进行 建 模 ,以 便 开 


MFR 各 种 工作 模式 的 运行 时 长 不 同 ,导致 截获 的 信 
号 样本 集 存在 工作 模式 样本 不 平衡 的 情况 ,也 给 
MFR 工作 模式 识别 增加 了 一 定 难 度 。 

综 上 ,针对 无 人 机 集群 截获 的 信号 样本 难以 直接 
融合 分 析 以 及 在 训练 样本 较 少 且 工作 模式 样本 不 平 
衡 时 难以 训练 出 有 效 分 类 器 的 问题 ,提出 了 一 种 基于 
Á HRES (self-organizing map ,简称 SOM) RE 
滑 图 信号 生成 的 MER 工作 模式 识别 方法 。 利 用 分 
布 式 SOM 算法 对 截获 的 信和 号 样本 进行 聚 类 ,使 得 同 


展 算 法 研究 。 如 表 1 所 示 , 一 个 雷达 短语 B 包含 4 
个 雷达 字 

B= [ww, wwal a0,cd E [1,9], (1) 
一 个 雷达 字 w HET m 

w= [pisos Pis" PL] ISLA L, (2) 
由 于 非 合作 方 通常 只 能 通过 分 析 MFR 信和 号 特征 来 
推测 其 意图 ,而 载 频 (RF) pki PW) .脉冲 重复 间隔 
(PRI) 能 提供 许多 有 关 MFR 状态 的 关键 信 
息 0 "中 ,将 脉冲 表述 为 
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I IV 口 ^H 
R1 “水 星 ”MFR 工作 模式 及 其 对 应 的 雷达 短语 p =L Ppr Prw Prde (3) 
工作 模式 雷达 短语 工作 模式 。 ”雷达 短语 其 中 , prr 、Ppw、Prr 分别 为 RF、.PW、PRI。 此 时 雷达 
Wiwawiws wiwrwwy 短语 可 重 写 为 
WoWiWswi W2W;W;W. B= [phopty Pmi], &—12,3444 (D 
ELLA WWW; Ww, 其 中 , phe 为 第 & 个 雷达 字 中 所 有 脉冲 的 RF 序列 。 
搜索 Wswiwaw Wawr ww 由 于 引入 无 人 机 集群 截获 MFR 信号 ,截获 的 全 部 信 
Wiwawaw WW; ww 号 样本 可 表示 为 
WsWsWiWs WeW7W7W7 X={B” | n=1,2,, N, m=1,2,.…,M,)}, 
Wawiwaw; Wiwswaws (5) 
www, Waws wsws 其 中 , B7 为 无 人 机 接收 的 第 m 个 雷达 短语 。 
Wawawawa Wawswsws 最 后 ,设计 一 种 MFR 工作 模式 识别 方法 fO). 
捕获 W3W3W3W3 W4WgWgWs 即 
Va 跟踪 保持 wswswsws Y=f(X), (6) 
A. d Ep, Y-(0"7[n21,2,7 Nas m-—1,2,75M,), 
v WiWgWgWg Wi1Wg9Wg9Wyg i — v ` 
Ee BU MM $t 为 无 人 机 7 接收 的 第 m 个 雷达 短语 的 分 类 结果 。 
A "i Yawswows 2 基于 SOM 聚 类 平滑 图 信号 生成 的 MFR 
a em WgWgWgWi W4WgWgWg 工作 模式 识别 方法 
le WowWee Ws WaoWowy - XE ER. 
N PUO MN 系统 框图 如 图 2 所 示 。 首 先 ,利用 无 人 机 集群 截 
C "v NAR Ik MFR 的 信号 样本 ,并 采用 分 布 式 SOM 算法 对 其 
C e MA HETRE EGG] — RAEE A I e REA RT REUR F 
O wa ib FI, ELFR A 63.6 (ci EH FE ui FUIT — s DL 
-elL 性 ,而 距离 较 远 的 簇 包含 的 信号 样本 差异 通常 较 大 。 
无 人 侦察 机 1 > a 
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1 无 人 机 集群 与 多 MFR 对 抗 场景 


第 1 期 


其 次 ,依据 上 述 聚 类 结果 的 特性 生成 关于 信号 样本 集 
的 平滑 图 信号 , 即 对 于 任 一 聚 类 复 内 的 样本 而 言 ,在 
图 信号 上 表现 为 禾 内 样本 全 连接 .与 相 邻 复 内 的 部 分 
样本 连接 ,使 得 图 中 任 一 节点 的 一 阶 邻 域内 通常 没有 
与 之 差异 较 大 的 节点 ,以 便 进 行 节 点 融合 处 理 。 最 
后 ,将 图 结构 以 及 各 图 节点 对 应 的 信号 样本 特征 作为 
GAT 的 输入 ,进行 节点 数据 融合 与 分 类 ,实现 MFR 
工作 模式 识别 。 在 数据 融合 阶段 仅 考虑 一 阶 邻 域 ,以 
尽 可 能 地 避免 节点 受到 非 同类 节点 的 影响 ,从 而 提取 
出 各 工作 模式 的 差异 性 ,应 对 工作 模式 样本 不 平衡 ; 


在 进行 数据 融合 的 同时 完成 类 别 信息 传递 ,应 对 训练 
样本 较 少 情况 。 


信号 特征 


2 ”系统 框图 


1; 基于 SOM 的 平滑 图 信号 生成 算法 


GAN 与 一 般 聚 类 算法 不 同 ,SOM 算法 依据 一 定 拓扑 

结 涛 初始 化 聚 类 簇 中 心 , 图 3 为 聚 类 篮 拓 扑 结构 ,图 
中 贿 贺 表示 到 类 秘 。 在 迄 代 过 程 中 , 先 计算 输入 信号 
桂 而 与 各 到 类 簇 中 心 的 距离 ,与 之 距离 最 小 的 簇 为 优 
肚 绎 ;然后 将 优胜 簇 及 其 邻居 簇 的 中 心 朝 该 信号 样本 
方向 调整 ,而 拓扑 结构 上 与 优胜 艇 距离 较 远 的 簇 基本 
不 受 影响 ,进而 保持 复 折 扑 结构 不 变 , 且 相 邻 复 包 含 
的 信号 样本 具有 较 强 的 相似 性 。 


图 3 聚 类 签 的 矩形 网 格 拓扑 结构 


采用 的 SOM 算法 适用 于 无 人 机 集群 分 布 式 处 
理 架 构 ,如 算法 1 所 示 。 Nw 为 以 第 4 行 第 v PURSE 
为 中 心 的 邻 域 函数 ， 
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算法 1 分 布 式 SOM 算法 

IMA: X = {B} | n=1,2,…,N,, m= 二 1,2,"…， 
M,) . 抢 形 网 格 的 行 数 r 和 列 数 c, 初始 学 习 率 qs 
初始 邻 域 因子 co, RRRA E. 

d Lee EAE Hu QOO — (4,00 | i =1,2, 
esr, jo—1,2,27550)5; REH S —(s" [n—1,2, 
es NS m —].2,*.M,). 

随机 初始 化 QO) 

t 40 

9) -— Mo 

olt) <o, 

for iter «- 1 to E do 

forn <— 1 to N, do 
form <— 1 to M, do 
[u,v] argmin(Bz — gq;,; Q0)) 
对 所 有 艇 执行 q GT D 4,0) H 
ONLA) [B7 — q; (1)] 


92 c D-—w/t 
oC +1) «o, /t 
pL. 
end for 
end for 


end for 
forn <— 1 to N, do 
form <— 1 to M, do 
[uv] --argmin(B — qi;,; 02) 
sm sce] 
end for 
end for 

对 于 SOM 聚 类 结果 而 言 ,拓扑 上 相 邻 的 聚 类 簇 
在 数据 空间 中 也 相 邻 , 即 拓扑 上 相 邻 的 复 所 包含 的 样 
本 通常 较为 相似 ,因此 依据 簇 邻居 关系 能 够 生成 关于 
言 号 样本 集 的 平滑 图 信号 。 根 据 簇 内 样本 全 连接 、 相 
邻 复 的 部 分 样本 连接 .不 相 邻 的 簇 不 作 处 理 的 原则 ， 
给 出 一 种 基于 SOM 的 平滑 图 信号 生成 算法 。 算 法 
WF: 

D 初始 化 邻接 矩阵 A — Lu, Es, WREN Ny 
HIRERE, Nx 为 样本 总 数 ,初始 化 邻居 簇 的 连接 
比例 为 re; 

2) 选择 一 个 复 为 主 簇 , 令 主 簇 所 含 样本 集 的 索 
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引 为 d, € R^". 其 中 Na 为 主 篮 内 样本 数 。 将 这 些 
样本 作 全 连接 处 理 , 即 邻接 矩阵 中 对 应 项 置 1， 
A(d,.d,)—1; 

3) 49 EIC 48 d Br FEAS 53 3k ffe op 
离 , 取 其 中 距离 最 小 的 若干 个 样本 的 索引 d, = 
[4 ,dd ,其 中 di € RN, 为 邻居 簇 
数 , ”为 第 ; 个 邻居 簇 所 含 样本 数 ; 

4) 从 d, 中 选取 一 个 值 , 令 其 为 ius 在 主 簇 中 随 
机 取 N n Xr 个 样本 的 索引 d 和 ,执行 A(i,,d') = 
l. 重复 该 步骤 ,直至 d, 中 所 有 项 都 执行 过 步骤 4); 

5) 返回 步骤 2» ,直至 所 有 簇 都 曾 被 设 为 主 徐 。 


2.2 图 注意 力 网 络 


OAT 由 若干 个 图 注意 力 层 堆 秋 构造 0 中。 任 一 
层 的 输入 为 邻接 矩阵 和 图 信号 中 所 有 节点 的 特征 向 
£d — [hh hy]; € RF PN 为 节点 总 
CST 为 节点 特征 维度 ,输入 层 的 图 节点 特征 向 量 如 
式 @) 所 示 。 为 实现 特征 提取 ,引入 一 个 可 学 习 的 线 
性 变换 矩阵 W ER ,此 时 节点 了 对 节点 ; 的 注意 
PRHE H 

e; = Attention (Wh; .Wh;) . (8) 
, AttentionC«) HJEM, Attention: R” 
RSR. 只 考虑 节点 一 阶 邻 居 对 节点 的 影响 ,并 
Fle; 作 softmax 归 一 化 ,得 


exp(e ;;) 
5 exp(ea) 


k kEN’ 
HIR: a; 为 节点 7 对 节点 i 的 注意 力 系 数 ; N 为 节 
点 QD 的 一 阶 邻居 节点 集 。 采 用 LeakyReLU 函数 增强 
非 线 性 ,并 将 Attention(Wh; ,Wh;) 展开 ,得 
exp(LeakyReLU(a* [Wh, || Wh,] )) 
5 exp(LeakyReLU(a" (Wh; || Wh,] )) i 


k€N' 


(9) 


a;; —softmax; (ej) = 


Qij 


a0) 
其 中 a € R2 Jj Attention) 的 权 值 向 量 。 此 时 单 
层 输出 为 
h'; =0( P 3 a;Mh;) ， (11) 
jEN’ 

其 中 , oC) 为 激活 函数 , h'; € RF 。 引 入 多 头 注意 

力 机 制 后 , 即 多 个 权重 , 单 层 的 输出 为 

K 
h';— | sc(2]a5Wh,). (42) 
Ee] 6" 

其 中 , K 为 注意 力 数 , 1 € RV, GAT 输出 层 通 
过 对 多 头 注意 力 取 平 均值 ,得 到 最 终结 果 "*1 , 即 信和 号 
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e HITTaA TVES T FRATI 


样本 所 属 的 工作 模式 ， 
" =o (> Safwn). a3) 


3 仿真 实验 与 分 析 


3.1 实验 设置 


仿真 场景 设 MER 数 N, = 3、 无 人 机 数 N, = 
100. 为 体现 一 般 性 ,将 无 人 机 随机 分 布 于 MFR 的 
扫描 空域 ,无 人 机 处 于 任 一 MFR 波束 内 即 认 为 该 
言 号 接收 成 功 。 每 次 仿真 截获 信号 样本 数 约 为 
3 500 个 ,不 平衡 度 ( 样 本 数 最 多 与 最 少 的 类 别 的 平 
均 样 本 数 比例 ) 约 为 10:1 ,假设 各 MFR 雷达 字 参 数 
如 表 2 所 示 。 为 降低 实验 结果 的 偶然 性 偏差 ,规定 
在 每 次 仿真 中 各 MER 均 从 最 大 参数 范围 内 随机 选 
定 特征 参数 的 子 区 间 ,特征 参数 为 RF, PW 和 PRI. 
并 将 该 子 区 间作 为 该 次 仿真 参数 范围 。 若 无 特别 
说 明 ,随机 变量 均 服 从 均匀 分 布 。 根 据 文献 L22] 设 
定 表 2 中 的 几 种 脉 间 调制 方式 ,其 中 国定 调制 方式 
表示 在 一 个 雷达 字 内 ,脉冲 特征 参数 值 固定 不 变 ， 
该 值 从 一 定 范围 内 随机 选取 ; 滑 变 调制 方式 下 , 雷 
达 字 中 各 脉冲 的 特征 参数 根据 脉冲 数目 等 间隔 地 
从 最 小 值 取 到 最 大 值 ; 捷 变 调制 方式 下 ,脉冲 特征 
参数 在 一 定 范围 内 随机 取 值 ;抖动 调制 方式 下 , 脉 
冲 特征 参数 围绕 某 一 中 心 值 在 十 5% 范 围 内 随机 拌 
动 。 另 外 ,为 使 各 雷达 字 时 长 基本 相同 , 令 雷 达 字 
w; ws 均 包 含 30 个 脉冲 ,wi 一 ws HEE 40 个 肪 
冲 ,w; 一 ws HEE 50 个 脉冲 。 


表 2 MER 雷达 字 参 数 设 置 


参数 WI Ws Wa wse W7 Wy 
最 大 参数 范围 /GHz [3,4] [4. 5] [5. 6] 
RF 子 区 间 长 度 /MHz 200 200 200 
脉 间 调制 方式 固定 滑 变 捷 变 
最 大 参数 范围 /ps [1. 2] [2. 3] [3, 4] 
PW  ” 子 区 间 长 度 /js 0.3 0.3 
脉 间 调制 方式 固定 抖动 滑 变 
最 大 参数 范围 /us [150.200] [100.150]  [50. 100] 
PRI ” 子 区 间 长 度 /us 10 10 
脉 间 调制 方式 固定 滑 变 抖动 


为 模拟 实际 环境 ,仿真 考虑 了 3 类 样本 缺陷 ,如 
表 3 所 示 。 由 于 PRI 是 相 邻 两 脉冲 到 达 时 间 (Time 
of Arrival. 简称 TOA) 的 差分 值 ,通过 对 TOA 添加 
测量 噪声 以 影响 PRI。 测 量 噪声 服从 均值 为 0, 标 准 
差 如 表 3 所 示 的 高 斯 分 布 。 


Y — dE] 


A3 不 同 场景 参数 设置 


RF,PW,TOA 
场景 ”测量 噪声 标准 差 / ”脉冲 丢失 率 /% ”干扰 脉冲 率 /% 
(MHz. ps. ps) 
1 [0. 0. 0] 0 0 
2 [2. 0.05, 2] 5 5 
3 [4, 0. 1, 4] 10 10 
4 [6. 0. 15. 6] 15 15 


SOM 算法 的 单行 . 列 聚 类 簇 个 数 均 为 30, 初 始 
学 习 率 y。 和 初始 领域 因子 oo 分 别 为 0.7 811.5. 4B 
居 艇 的 连接 比例 x. 为 0.01, 邻 居 簇 的 数目 为 4, 即 
E FÆ ATHERE. KRH GAT 包 
含 2 个 图 注意 力 层 ,第 一 层 的 注意 力 头 数 为 8, 每 头 
注意 力 提取 的 特征 数 为 20, 然 后 利用 ELU(Cexponen- 
tial-linear unit) 激 活 孔 数 增强 非 线 性 ;第 二 层 用 于 分 
类 = 了 注意 力 头 数 为 1, 激 活 函 数 采 用 softmax, 注意 力 

的 特征 数 为 15 ,分别 对 应 15 个 类 别 , 学 习 方式 
为 喜 推 式 学 习 , 和 迭代 轮 次 为 1 000, 优 化 器 为 Adam, 

意 力 层 的 dropout 比例 为 0. 2。 本 算法 在 计算 距 
离 时 采用 欧 氏 距 离 。 为 消除 量 纲 影响 ,对 数据 进行 标 


一 e 一 未 方法 的 准确 率 
Fa% —e— 本 方法 的 Fl 指数 0.2 

m 一 一 HSSLSTM 的 准确 率 || 0 ] 
Qo. —4— HSSLSTM 的 F1 指 数 || O 

012345 10 15 20 25 
每 种 类 别 训练 样本 数 
(a) 场景 1 

1.0 1.0 

0.9 0.9 

0.8 10.8 

0.7 0.7 

0.6 10.6 
€ 0.5 0.5 3E 
E: E 

0.4 104 

0.3 103 

一 @ 一 本 方法 的 准确 率 
0.2 一 e 一 本 方法 的 F1 指 数 0.2 
一 和 -HSSLSTM 的 准确 率 
0.1 一 HSSLSTM 的 Fl 指数 |] 0.1 
012345 Tio 15 20 25 
每 种 类 别 训练 样本 数 
(0 场景 3 


F SOM 聚 类 平滑 图 信号 生成 的 MER 工作 模式 识别 方法 


准确 率 


准 化 处 理 . 
7 mean(Cx) 
stdCz) 
其 中 , meane) std 分 别 为 均值 标准 差 。 硬 件 配 
置 为 :Intel@ Core™i9-10900X CPU (22. 20 GHz 中 
KAFEE; NVIDIA GeForce RTX 3070 图 像 处理 器 ; 
16 GiB 内 存 。 


3.2 方法 性 能 对 比分 析 


为 验证 本 文 方法 的 性 能 ,对 比 了 基于 分 层 序列 到 
序列 长 短期 记忆 网 络 的 MFR 工作 模式 识别 方法 ( 记 
为 HSSLSTM) , 除 特别 说 明 外 ,该 算法 的 仿真 参数 
与 文献 [15] 一 致 。 为 讨论 在 训练 样本 较 少 时 的 识别 
效果 ,仿真 了 每 种 类 别 训练 样本 数 m. — 1.2.3. 4.5. 
10,15,20,25 的 情况 ;为 讨论 样本 缺陷 对 识别 结果 的 
影响 ,分 别 对 表 3 中 4 种 场景 进行 实验 。 采 用 准确 
率 、F1l 指数 (常用 于 类 别 不 平衡 场景 ) 作 为 评价 指标 ， 


, (14) 


上 述 每 种 情况 都 进行 100 次 仿真 实验 并 取 其 结果 平 
均值 ,如 图 4 所 示 。 
xx 
y 
E 
一 e 一 本 方法 的 准确 率 
一 @ 一 本 方法 的 F1 指 数 
一 上 一 HSSLSTM 的 准确 率 
一 9 一 HSSLSTM 的 F1 指 数 
012345 10 15 20 25 
每 种 类 别 训练 样本 数 
(b) 场景 2 
0.9 0.9 
0.8 0.8 
0.7 10.7 
0.6 10.6 
0.5 10.5 $ 
yz 
0.4 Q4 E 


0.3 103 
一 @ 一 本 方法 的 准确 率 
一 9 一 本 方法 的 Fl 指数 02 
0.1 一 一 HSSLSTM 的 准确 率 | | 0 1 
一 9 一 HSSLSTM 的 Fl 指数 | ] 
012345 10 m — ; 
每 种 类 别 训练 样本 数 
(d) 场景 4 


4 2 种 方法 在 4 种 场景 下 的 识别 结果 
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由 图 4 可 知 , 随 着 训练 样本 数 的 增多 ,各 场景 下 
2 种 方法 的 准确 率 和 FI 指数 均 呈 增长 趋势 。 由 于 
HSSLSTM 采用 的 神经 网 络 模型 相对 较为 复杂 , 需 
要 优化 的 参数 较 多 ,因此 在 训练 样本 较 少时 泛 化 性 能 
较 差 ,尤其 在 n, < 5 时 ,基本 无 法 训练 出 有 效 分 类 
器 。 而 本 文 方法 通过 生成 关于 信号 样本 集 的 平滑 图 
信号 ,将 相同 工作 模式 下 的 信号 样本 进行 连接 ;然后 
依据 样本 之 间 的 关联 ,采用 GAT 模型 进行 节点 信息 
融合 与 分 类 。 其 中 ,将 融合 范围 限制 在 一 阶 邻 域 ,从 
而 减少 不 同 工 作 模 式样 本 间 的 干扰 ,以 应 对 工作 模式 
样本 不 平衡 ;利用 GAT 进行 邻 域 类 别 信息 传递 , 以 
摆脱 对 大 规模 训练 样本 的 依赖 。 

对 于 场景 1 , 即 图 4(a) 所 示 。 当 72. 一 25 时 , 即 训 
练 样本 数 占 比 约 为 10.71%, 本 方法 的 准确 率 和 Fl 
指数 分 别 为 98. 95%、98. 39% ,而 HSSLSTM 识别 准 
WERA Fl 指数 分 别 为 76. 1596. 76. 05% , 4H EC 43 30] 
BEDT 22.894,22. 34 外 。 但 本 方法 运行 时 间 相对 于 
HSSLSTM AK n, —25 时 ,本 方法 中 生成 平滑 图 
信号 的 平均 用 时 为 2. 11 s, 训 练 GAT 的 平均 用 时 为 
29413 s. 而 HSSLSTM 的 平均 用 时 为 7.52 s. 
HSSLSTM 用 时 较 少 的 主要 原因 是 训练 样本 数 少 ， 

LGAT 需要 进行 节点 融合 操作 , 故 用 时 相对 较 多 。 
四 通过 对 比 4 种 场景 下 2 种 方法 的 识别 准确 率 和 
险 糙 数 可 发 现 , 随 着 场景 中 噪声 干扰 愈 发 严重 ,2 种 
方 潜 的 识别 准确 率 和 F1 指数 基本 呈 下 降 趋势 。 这 
已 遇 于 截获 的 信号 样本 受 噪声 干扰 导致 存在 一 定 程 
度 久 样本 缺陷 ,利用 这 些 存在 缺陷 的 样本 训练 算法 模 
型 稳 降 低 模 型 的 泛 化 性 能 。 但 即使 在 场景 4 中 , 即 图 
4( 册 所 示 , 本 方法 在 .一 25 时 的 识别 准确 率 和 Fl 指 
数 分 别 为 84. 3896.82. 19% ,说明 本 方法 能 够 克服 一 
定 程度 的 样本 缺陷 。 

3.3， 聚 类 数 对 识别 结果 的 影响 分 析 


为 了 探究 分 布 式 SOM 算法 的 聚 类 数 对 MFR T. 
作 模 式 识别 效果 的 影响 ,仿真 了 6 种 单行 . 列 聚 类 数 
取 值 下 本 方法 的 识别 准确 率 ,每 种 类 别 训练 样本 数 为 
25。 此 外 ,还 考虑 了 表 3 所 示 的 4 种 场景 ,对 每 种 情 
况 进 行 100 次 仿真 实验 并 取 结 果 平 均值 ,仿真 结果 如 
表 4 所 示 。 

从 表 4 可 看 出 ,识别 准确 率 随 聚 类 数 增多 而 提 
高 ,但 在 不 同 噪声 干扰 条 件 下 影响 不 同 ,在 噪声 干扰 
较 小 时 ,如 场景 1 和 场景 2, 单 行列 聚 类 复数 在 大 于 
15 时 的 增长 率 开 始 趋 于 平缓 ; 而 在 噪声 干扰 较 大 时 ， 
如 场景 3 和 场景 4, 这 一 增长 率 依然 较 高 。 具 体 而 
言 ,在 场景 1 下 ,单行 、. 列 聚 类 复数 为 15 和 30 时 的 识 


| .,,..2023 年 2 月 
T IAT FEATH 
表 4 NARAT MFR 工作 模式 识别 准确 率 


L aM 2 
LIAA IV 


单行 、 列 复数 ”场景 1 场景 2 场景 3 场景 4 
5 0. 888 1 0. 814 6 0. 7011 0. 603 2 

10 0. 950 3 0. 8913 0. 774 5 0.695 8 

15 0. 977 3 0.919 8 0.8214 0. 740 7 

20 0. 986 6 0. 944 3 0.872 5 0.812 7 

25 0. 987 6 0. 957 2 0. 898 3 0. 829 7 

30 0. 989 5 0. 968 9 0. 916 3 0.843 8 


别 准确 率 仅 相 差 1. 2276s 而 在 场景 4 下, 单行, 列 聚 


类 复数 为 15 和 30 时 的 识别 准确 率 相差 10. 3176, 
这 是 由 于 场景 1 中 未 考虑 噪声 干扰 ,因此 依靠 少量 聚 
类 艇 即 可 生成 较 好 的 平滑 图 信号 , 且 GAT 的 分 类 性 
能 也 未 受到 噪声 影响 ; 而 在 场景 4 下 ,由 于 噪声 干扰 
引起 的 样本 缺陷 较为 严重 ,使 得 生成 平滑 图 信和 号 
GAT 训练 过 程 都 受到 较 大 影响 ,导致 识别 精度 显著 
下 降 。 


4 结束 语 


从 建立 无 人 机 集群 截获 的 MFR 信和 号 样本 的 融 
合 依据 出 发 ,针对 现 有 方法 在 训练 样本 较 少 且 工作 模 
式样 本 不 平衡 情况 下 识别 精度 低 的 问题 ,提出 了 一 种 
基于 SOM 聚 类 平滑 图 信号 生成 的 MFR 工作 模式 识 
别 方法 。 利 用 分 布 式 SOM 算法 对 无 人 机 集群 截获 
的 信号 样本 集 进 行 聚 类 ,依据 聚 类 结果 以 平滑 图 信和 号 
的 方式 表征 信号 样本 集 , 采 用 GAT 完成 图 节点 分 
类 ,实现 MFR 工作 模式 识别 。 实 验 结果 证 明了 本 方 
法 的 可 行 性 及 对 于 样本 缺陷 问题 的 适用 性 。 但 本 方 
法 在 应 用 分 类 器 阶段 未 考虑 分 布 式 架 构 , 针 对 该 问题 
还 需 进 一 步 研 究 。 
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